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1. Introduction 

DifFractive deep inelastic scattering, or "hard diffraction", is a subclass of semi- inclusive deep 
inelastic reactions in which a virtual point-like particle (quark or lepton) is assumed to interact 
directly with the pomeron emitted by the proton target. This may happen in a special kinematical 
configuration when the hit proton continues its motion nearly in the forward direction. 

There are several reasons why this class of reactions is receiving so much attention. The first 
and most obvious one comes from the prospect to study the internal structure of the pomeron, 
a hypothetical quasi-particle with the quantum numbers of vacuum, responsible for diffraction in 
the Regge pole model. Although the composite nature of the pomeron, within the context of the 
quark model and QCD, was never a question for the theorists, the possibility of its experimental 
verification in deep inelastic scattering was first formulated in Ref. j^] . In subsequent experiments 
of the UA8 Collaboration 2-jet events were reported as evidence of the partonic structure of 
the pomeron. 

The study of the so-called "large rapidity gap" events at HERA initiated by the ZEUS Col- 
laboration [||, revitalized the subject. HERA produced a large number of high-precision data 
in a wide kinematical range of X and published recently by the HI g and ZEUS [| Col- 
laborations. A new and important development in the subject is the experimental study of the t 
dependence of the cross-sections, ignored until recently. 

Theoretical studies of the diffractive deep inelastic scattering (hard diffraction) were pioneered 
by the papers of Donnachie and Landshoff long before the HERA experimentations. Now there 
is a large literature (partially listed in Refs. ||7|"[l3||, we apologize to those not included because 
of the limited scope of the present paper) dealing with various aspects of the phenomenon. In our 
opinion, the present state of the subject can be summarized as follows. 

There is nothing surprising in the very existence of the large rapidity gap events (occurring e.g. 
also in the pp diffraction dissociation) provided diffraction, the Regge pole model and properties 
of the S matrix do not change essentially when one of the external particles goes off shell. The 



last assumption is the most delicate one in the whole subject, and we shell come back to it. 

The next question is whether the appearance of a rapidity gap, within which secondary particles 
are not produced, is a manifestation of diffraction, and if so: is it a pomeron exchange, or can 
it be simulated by an alternative mechanism ?. The first answer is "yes" - to the extent other 
features of diffraction, namely typical <^ (or a;p) and t dependence will be confirmed. 

The ^ dependence has been recently measured and found to be « with the values a — 
1.30 ± OM{stat) """H'^^ {sys) as measured by ZEUS |] and a = 1.19 ± 0.06{stat) ± 0.07{sys) by 
HI These two values are compatible and both agree with the pomeron intercept ||l^ most 
reliably extracted from the p — p and p — p total cross section. The t dependence is now being 
measured and we discuss it below. 

Whether diffraction is mediated by a pomeron (pole) or not is a matter of convention. A 
factorizable vacuum Regge pole exchange with the trajectory whose intercept is slightly beyond 1, 
a(0) = 1 + e, is known to be the adequate mechanism of diffraction, at least in hadronic reactions. 
Whether and how is the formalism affected by off mass shell effects is an open question. In any case 
the use of the pomeron exchange is more justified in diffractive DIS, where at least the pomeron 
interaction is directly related to hadronic diffraction, than generally in a (small-x) DIS, related to 
the absorptive part of the elastic "f* — p scattering only for < = 0. 

The large flow along the photon line was an argument for many authors to use perturbative 
QCD calculations for a single- ^ or two-gluon (pomeron) exchange in diffractive DIS. In 
our opinion, the use of perturbative QCD may be justified only in resolving the pomeron struc- 
ture (QCD evaluation of the photon pomeron vertex), while the exchanged object, the pomeron 
propagator, is essentially nonperturbative and it is the same, whatever is the photon virtuality. 
A more detailed discussion concerning the uniqueness of the pomeron in elastic and deep inelastic 
scattering may be found in Ref. |p^ . 
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The above construction relies very much on the hypotheses of factorization of the diffractive 
structure function F^^'^'' into a product of the pomeron flux and the pomeron structure function. 
SmaU deviations from factorization, e.g. due to muhipomeron exchanges, initial and final state 
interactions etc. are admissible and certainly should be accounted for in future calculations. At 
the moment, according to the recent measurements at HERA |^,||, factorization is confirmed at 
the present level of accuracy. 

The ultimate goal of the experiments at HERA and of their theoretical interpretation is to an- 
swer the question: what is the pomeron structure or, in other words, to find the right parametriza- 
tion for the pomeron structure function. Because of the complexity of the problem and the existing 
uncertainties, both on the experimental and theoretical side, the final result is still ambiguous and 
sensitive to the input. 

With the present paper, we try to contribute by one more step in the clarification of the 
remaining uncertainties. First of all in Sec. 2 we introduce notations and define the kinematics. 
Next in Sec. 3 we discuss the pomeron fiux and its (re)normalization. The pomeron structure 
function and the comparison with the experimental data are developed in Sec. 4 while Sec. 5 is 
devoted to a study of the t dependence of the distribution function. We use the earlier experience 
based on the description of hadron diffraction with a non-linear pomeron trajectory to predict 
the t dependence of the diffractive DIS. 

2. Notations and kinematics. 

The notations and kinematics for the process 

e-(fc) +pip) e-(fc') +pip')+Xipx) (1) 

are standard. The fourmomenta of the virtual photon and of the exchanged reggeon are q — k~k' 
and r = p — p' , respectively. Besides the usual DIS variables = —q^^x — Q^/{2p ■ q), and 
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— {p + q)'^ the new variables 



r ■ q X 



2r • p ■ q j3 

are introduced. Sometimes ^ is called also xp. 



(2) 



Let Mx be the invariant mass of the hadronic system X, = ^'^x^ then the fourmomentum 
transfer squared, t — , ranges between the extreme limits 



2W 



2W 



1/2 



2VK 



M 



X 



1/2' 



(3) 



or 



t_ ~ — mpO;^ ( 1 



X 



M_ 



2 /-2 



while t+ ~ —W'^ is a formal limit since \t\ is less than ^ 7 GeV'^ for the measurement in 
Experimentally, the i distribution has never been measured until recently because of the existing 
difficulties in identifying the proton hit by the photon but continuing its motion in the nearly 
forward direction. 

On the other hand it is well known from the hadronic physics that diffraction is typical of the 
domain of about \t\ < 2GeV^. In this region, the elastic differential cross section is known to 
decrease almost exponentially up to about \t\ ~ 1 GeV^ , followed by a dip-bump structure between 
1 and 2 GeV^ . The latter is an important feature of high energy diffraction. The deviation from 
the exponential behaviour of the cone may be taken into account by a nonlinear trajectory, a 
simple representative example of which is jl^ 



a{t) — aQ + ait — a2 ln(l — a^t) 



(4) 



where = 1 + e and ai(i = 1, 2, 3) are parameters. 



The diffractive structure function is defined from the cross section as in [g| and the factorization 
hypothesis can be written as 

F]P/pi£,it) and G^^jp{f3,Q^) are, respectively, the pomeron flux and the pomeron structure func- 
tion, to be introduced in the next two sections. We quote also the upper and lower kincmatical 
limits for ^ typical of the HI and ZEUS data: 





HI 


ZEUS 




0.05 


0.01 




3 X lO-'^ 


6.3 X 10-4 



3. The pomeron flux and its (re)normalization 

As in Ref. we fix the form of the scattering amplitude following the duality prescription. 

In dual models the dependence on the Mandelstam variable t enters the amplitude only through 
the trajectory a(t). Hence in the Regge limit, s — > cx) at fixed t, 

^(s,i) =e^(^)"(*) , (6) 

where B{s) = B^i + ln(s/so) ^ i7r/2. 

We have already discussed the relevance of a nonlinear trajectory that extrapolates between 
"soft" and "hard" scattering. The interest for a possible flattening of the pomeron trajectory 
increased in view of the recent experimental result of the UA8 Collaboration In order to 

make clear the effect of a nonlinear a{t) we will consider Eq. (^) in the two limiting cases 

a) a2 = : = a(0) + a'i , (7) 
b) ai = 0, a2 = 0^3 : Oi{t) — a{0) — 7 ln(l — -ft) . (8) 
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The first instance has been studied in Ref. for different choices of a{0) and a' — 0.25. In 

case b) we will choose 7 0.5 in order to reproduce the same slope near t = and to deal with a 
trajectory close to the optimal solution of Ref. 

The constant term in -B(s), B^i, has been determined from a fit to the clastic p — p differential 
cross section at ISR . For the pomeron flux 

Fp/p = Ce2(^-'"«)"W.e (9) 

we get 

B — Bdiff = Bel/2 2± 7.0. 

The main hypothesis here is that only the pomeron trajectory contributes to the single diffractive 
cross section, a condition supported from the experimental evidence for factorization. Since the 
trajectory (H) is an effective one, the best determination of a{0) comes from the overall Regge fit 
of total cross sections in Ref. [Q 

a(0) ~ 1.08. 

We will use in the following this value that differs slightly from the CDF intercept a{0) = 1.112 ± 
0.013 obtained from total cross section data, but provides a sensible determination of the 
constant in front of the flux (H). 

From [1^ ] we get the pomeron contribution to p — p total cross section and find 

c = Ce2S«(°) = —m5 = 1.1087 (10) 
167r 

From the same fit we deduce that the contribution to the cross section of non-asymptotic Regge 
terms {p, /,...) is rather small at HERA, approximately 8% at the lowest W values. As seen 
in Fig. 1, there is a marked difference in the t dependence of the pomeron flux derived from a 
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logarithmic trajectory, curve (1), or a linear one, curve (2). Curves (3), ||^ and (4), ^ refer to 
models for the pomeron flux approaching the limiting cases represented by eqs. (Q) and (|^). 



The integrated pomeron flux for the nonlinear trajectory ^ is 



c (l~7^_)-^(g)-(l-7M-''^^\ .l-2o(0) ^ ^.1-2.(0) .1 . ^ 



where 



&(0 = 27(^3 -In 0-1. 

The approximation in ( pT| ) is reasonable since t_ cx — and \t+\ is rather large, ~ 7 GeV'^ 
for the measurement described in Ref. Q. 

Now the ^, t dependence of the diffractive structure function is completely fixed within the 
model. In order to set the normalization and to make a comparison with the experimental data 
we need the pomeron structure function. 

Renormalization of the pomeron flux [T^ affects only the determination of the pomeron struc- 
ture function. If we require that no more than one pomeron should be exchanged in one diffractive 
proton interaction, then a bound must be imposed on the integral 

j^" ^ ^ e-^(i){£;z[/(a)] - E^[fi^H)]} , (12) 



where 



/(0 = 2(a(0)-l)(i?-;l-lnO 
27 

and S^H, £,L are the upper and lower kinematical limits for ^ specified in Sec. 2. In Eq. (|l2|), 
Ei(z) is the exponential integral and the numerical values of the integrated flux are 1.34 



for ZEUS and 2.42 for HI For a linear trajectory (|^) these values are somewhat smaller: 
1.23 and 2.23 respectively. Considering the unitarization procedure explained in Ref. flux 
renormalization should be applied also to some {Q'^,(3) bins of the ZEUS data, in particular for 
small and large f3. Since corrections remain within experimental errors we do not scale down 
the integrated flux. 

4. The pomeron structure function and comparison with experimental data. 

As in our previous paper we consider a pomeron composed mainly of gluons. This point of 
view reflects the structure of the BFKL pomeron ||2|] that, at least for large momentum transfers, 
relies on a sound basis. Non perturbative effects will represent a new contribution to the aforesaid 
picture but should not change the pomeron content. On its gluon content both theoretical [ p^ 
and experimental p3[ | constraints exist. Experimental data indicate that a large fraction of the 
pomeron momentum, up to 80%, can be carried by hard gluons [ p3[ and that the pomeron 
structure function is approximately independent of |^,|| . 

Since only quarks interact with the photon, quarks must be present in the input distribution 
with a fraction of the pomeron momentum larger than 20%. In Ref. we noticed that, at large 
/?, to the assumed gluon distribution 

pG{P,Q') = a{Q^){l~(3) (13) 

one must add the qq sea contribution. Quark loops and the mesonic trajectories, included in the 
effective trajectory (^), will give rise to a quark distribution 

Pq,{f3,Q')^d{Q')f3{l-/3) (14) 

for the quark i, whose form we borrow from Ref. Both distributions contribute at the starting 
scale Qq = 5 GeV^ and, due to the presence of quarks, the previous estimate of the evolution JTzt 
does not hold now. 



We can keep however the calculation simple and transparent if we neglect gluon recombination 
effects and limit ourselves to the low-Q^ region, < 40 GeV^. This condition will not destroy 
the predictive power of the result since "there is no evidence for any substantial dependence 
of F^" a function proportional to the pomeron structure function. 

We use a recursive method for solving the massless inhomogeneous Altarelli-Parisi equations. 
This method is based on the power expansion of the solution in the parameter s: 

and, since it has been explained at full length in Appendix A of Ref. ||2^, we will not repeat 
the details here. We choose Qq — 5GeV^,A — 0.2 GeV and, taking into account gluon to quark 
conversion, the result at small s is 

Pq,{P,Q^)^dPil-P)il-s) + 
+ es |^d/3(l + In (^^^) ) + ^« - + y + /Jln/S) } (15) 

with 

a = a{Ql), d = d{Ql) 

and 

6 

^ " 33-2/ 

where / is the number of quark flavors. 

With three flavors, f — ?>, the pomeron structure function is 

G,/p(/3,Q') = ^/3ft(/3,Q'). (16) 
Parameters a and d are not independent since we impose the momentum sum rule in the form 
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/' 



dp (3 

Jo 

obtaining the constraint 



d=l{2-a) (17) 



While a proof of the validity of this sum rule for the pomeron is lacking |25|, its applicability 



appears reasonable once a model for the pomeron in terms of its constituent is assumed |7y8|,|12| . 

Insisting on the particle nature of the pomeron we can evaluate the total momentum carried by 
quarks and antiquarks 

M{Q^)^6iyq,{/3,Q^)d(3 = d+^(^a-^-^d^ 

and, taking into account the relation (pT[), we find that the dependence disappear for a = 
154/95. 

This value is not far from the one obtained in fitting to the data, 

a = 1.458, 

that gives a weak dependence on with dM{Q^)/dQ^ < 0. The relative contribution of quarks 
to the pomeron momentum is near 0.25 in the range of where the recursive method applies. 
The presence of quarks at every scale simulates a quarkball |l3| . 

A plot of all ZEUS data for F^^^^ ^ for different /3 values, regardless their value, as in 
figures (2a-2c), shows that the dependence is indeed weak. The large errors permit only to 
say that our model is compatible with data. Continuous and dashed lines are the result of the 
evolution at = 16 and 28GeV^, respectively, with the nonlinear trajectory (||) in the flux. No 
visible change is noticed if the calculation is repeated with the linear trajectory (Q). 

In comparison the HI data, presented in Fig. (2d) for one j3 value show a larger spread in Q^. 
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In this case the selected values of for the fit are = 
(dashed curve) 



12GeV^ (continuous curve) and 25 6*6^^ 



5. Predictions and conclusions. 

In order to get a significant picture of the predictions from the model, we integrate over (3 and 
the differential cross section expressed as 



d(3dQ'^d£,dt 



[1 + (1 - vf] %/p(e, i)G^/p (/?, Q^) , (18) 



where 



and s — (296 061^)^. The /3 range of integration is restricted to the interval 

0.02 < /3 < 0.8 

and the upper limit avoids the region near /3 = 1 where theoretical and experimental uncertainties 
are larger. 

The integration has been performed for two different intervals, both comprised in the region 
less sensitive to theoretical approximations. The starting scale = 5 GeV^ for evolution is the 
lower limit in both cases and the d^a/d^dt results are displayed in Figs. 3a, 3b. Predictions based 
on the logarithmic trajectory (H), continuous lines in Fig. 3, are clearly distinguishable from the 
ones obtained with the linear trajectory (^), dashed lines. Future experimental data, even at 
moderate \t\ values, will be able to decide between the two possibilities. 

The t dependence shown in Figs. 1, 3a, 3b has a nature typically diffractive, i.e. it is strongly 
peaked in the forward direction. At present, the (nearly) exponential decrease of the differential 
cross section d^a/d^dt to large extent has been fed in by the choice of the residue and form of 
the pomeron trajectory. Nevertherless, this result is far from being trivial since in the formalism 
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under consideration, the t dependence is correlated with other variables and ultimately the choice 
of different inputs (pomeron trajectories) will be tested experimentally. 

It is well known - both from the S-matrix theory p6t and from the fits to hadronic data - 
that the pomeron trajectory contains a significant non-linear part. The amount and form of this 
non-linear correction is a matter of debate in the literature. Our choice of the pomeron trajectory 
(Q) was based on earlier fits ||l6| to the data on high energy elastic hadron scattering. Data 
on the t dependence of the diffractive structure function have been discussed recently by two 



experimental groups: ZEUS |^ and by UA8 |18|. We just notice that the UA8 results show 



the same trend in the behaviour of the differential cross section that derives from the use of the 
logarithmic trajectory (Q). The results of the measurement of the t dependence at HERA are 
preliminary and still not complete. 

The knowledge of the t distribution will settle important questions. For example, the photon- 
pomeron vertex may be also t dependent and it may bias the discussed behaviour of the structure 
function. 

Another prominent feature in the t dependence of diffraction, well known in elastic hadron 
scattering, is the apperance of the diffraction minimum. The expected effect should be visible 
"by eye" because of its unmistakable structure, but the energy (here, W) and t values where the 
minimum should appear are near the kinematical boundary of the relevant experiments on deep 
inelastic diffractive scattering. For the above (kinematical) reason it has not yet been seen even 
in hadronic diffractive dissociation. Our model, as well as others, does not contain this structure. 
Its experimental observation however would resolve all doubts concerning the diffractive nature of 
the "large rapidity gap events" and related phenomena. 
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Figure Captions 

Fig. 1: The pomeron flux -Fjp/p(Cj t) versus —t for ^ = 0.0032 with the logarithmic trajectory (1) 
or the linear one (2) discussed in the text. For comparison the pomeron fluxes of Refs. 0, 
(3) and (4) are also drawn. 

Fig. 2: The diflractive structure function F^^'^^P, Q'^ , x) versus x, for different values of P and 
g2 compared with the ZEUS data || (a), (b), (c) and with the HI data g, (d). F^^^'^ has 
been evolved at ~ 16 GeV^ (continuous line) and = 28 GeV^ (dashed line) for the 
ZEUS data. For the HI data, = 12GeV'^ (continuous hue) and = 25GeV'^ (dashed 
line). 

Fig. 3: The diflractive cross section (Pa/d£,dt predicted from the model, versus —t for diflerent ^ 
values. The order of the ^ values is respected in the curves. Two different integration regions 
for Q2 are considered: (a) 5GeV^ < < 20GeV^ and (b) 5GeV^ < < iOGeV^. 
Predictions from the logarithmic trajectory (full lines) and from the linear one (dashed 
lines) are shown. 
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